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PACS 42.65.Tg - Optical solitons; nonlinear guided waves 

PACS 42 . 65 . Jx - Beam trapping, self-focusing and defocusing; self-phase modulation 

PACS 42 . 65 . Wi - Nonlinear waveguides 

Abstract. - We report dissipative surface solitons forming at the interface between a semi-infinite 
lattice and a homogeneous Kerr medium. The solitons exist due to balance between amplification 
in the near-surface lattice channel and two-photon absorption. The stable dissipative surface 
solitons exist in both focusing and defocusing media, when propagation constants of corresponding 
states fall into a total semi-infinite and or into one of total finite gaps of the spectrum (i.e. in a 
domain where propagation of linear waves is inhibited for the both media). In a general situation, 
the surface solitons form when amplification coefficient exceeds threshold value. When a soliton 
is formed in a total finite gap there exists also the upper limit for the linear gain. 



Shallow periodic modulations of a refractive index may 
considerably affect the diffraction of low-power beams and 
properties of stationary nonlinear excitations [I] . Of par- 
ticular interest are truncated periodic structures that sup- 
port so-called surface lattice solitons propagating along 
the interfaces between periodic and uniform media. While 
the history of linear surface modes goes back to the semi- 
nal paper of Tamm [5] , one-dimensional surface modes in 
nonlinear lattices were predicted in [3] and subsequently 
obtained in both focusing [1] and defocusing [S1[B] mate- 
rials, while their two-dimensional counterparts were gen- 
erated in [J] . A characteristic feature of surface solitons 
is that they usually exist above the energy flow threshold. 
Surface lattice solitons inherit some properties of solitons 
in bulk lattices, in a sense that their shapes and trans- 
verse extent depend on the location of propagation con- 
stant in a band-gap lattice spectrum. Moreover, there ex- 
ists no sharp transition between modes strongly localized 
at the surface and modes localized at some distance from 
the surface or bulk modes: the localization domain moves 
smoothly from the boundary to some domain in the bulk 
as one follows the one-parametric branch of the solutions 
j8]. A surface can also support breather solutions, i.e. lo- 
calized modes with periodically changing shapes [8] . Even 
more rich behavior is observed when a defect is added into 



otherwise periodic guiding structure. In particular, when 
a defect is introduced into a surface lattice channel it may 
dramatically affect the conditions for surface soliton for- 
mation and result in appearance of new types of surface 
modes [5]. 

The above referred studies were dedicated to conser- 
vative systems. An interesting and still open problem is 
the impact of non-conservative surfaces on existence and 
properties of the localized modes. A particularly interest- 
ing situation arises when a gain is applied at the interface 
between a homogeneous and a periodic medium. It is rel- 
evant to mention that the effect of a localized gain was 
analyzed for existence and interaction of gap solitons in 
shallow fiber Bragg gratings [TU], in media of two-level 
atoms [llj , for formation of spatial dissipative solitons in 
systems governed by complex Ginzburg-Landau equation 
[12j . and for dynamic emission of moving lattice solitons 
in systems without dissipation [13] . 

In this Communication we address the properties of 
truly stationary dissipative solitons forming at the edge 
of semi-infinite lattice and existing due to the balance be- 
tween localized gain in the near-surface lattice channel and 
strong two-photon absorption in a cubic medium. Dissi- 
pative surface solitons [14] were studied before only in the 
truncated discrete systems governed by Ginzburg-Landau 
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equation with uniform gain (the model introduced in |15j). 
Here we show that in experimentally realistic setting dissi- 
pative surface solitons are pinned by the "defect" channel 
where gain is realized. Such modes are attractors and thus 
can be excited from sufficiently large class of initial condi- 
tions, ranging from localized or extended regular patterns 
to noisy patterns. We argue that surface solitons form 
when gain coefficient exceeds a threshold value, and if the 
propagation constant belongs to a finite gap, there exist 
also the upper bound for the allowed gain coefficient at 
which surface solitons can still be found. 

Specifically, we consider propagation of laser radia- 
tion in a lattice with spatially localized linear gain and 
nonlinear losses that can be described by the nonlinear 
Schrodinger (NLS) equation for the dimcnsionlcss light 
field amplitude q: 

m = -^m ~ i R M - '^( 7 ?)] i - cr l?l 2 '? _ ia k\ 2 Q- (i) 

Here 7/ = x/xo and £ = z/Ldif are the transverse and 
longitudinal coordinates normalized to the characteristic 
beam width xo and to the diffraction length Ldif = kx^, 
k = 27rno/A, A is the wavenumber, no is the unper- 
turbed refractive index, Iq is the characteristic inten- 
sity, L„i = no/kn2lo is the nonlinear self-action length, 
Pi = L^f I Lgainis linear gain coefficient, L ga i n is the am- 
plification length, a = Ldif/Li oss is the coefficient of non- 
linear losses, Li oss = I/012I0 characterizes the length of 
two-photon absorption and a = 1 (u = — 1) corresponds 
to focusing (defocusing) nonlincarity. 

Spatial solitons were successfully observed in nonlin- 
ear waveguide arrays made of AlGaAs alloy below the 
half-band-gap, A w 1.53 /im. Such arrays have typical 
length of 6 mm, waveguide separation of 4 -=- 7 /im and 
effective core area about 20 /im 2 [16]. In this case (as 
well as in our model) the two-photon absorption is the 
dominating mechanism of optical losses: the linear ab- 
sorption coefficient is around 0.1 cm -1 while the typi- 
cal value of two-photon absorption coefficient is ot2 = 0.3 
cm/GW, and the representative peak soliton intensity is 
above Iq « 5 GW/cm 2 . Importantly, the same semicon- 
ductor material (AlGaAs) is widely used for production 
of the wide-band (typical bandwidth of 60 -=- 70 nm) semi- 
conductor optical amplifiers with rather high optical gain 
in the same spectral range [17]. In particular, for the 
above mentioned structure, no = 3.34 and the Kerr co- 
efficient 17,2 = 1-6 • 10~ 13 cm 2 /W. Then for a laser beam 
with xq = 5 /im one obtains Ldif ~ L n i ra 0.34 mm and 
Lioss ~ 6.67 mm, that corresponds toa~ 0.05 (TE polar- 
ization is under consideration). Then pi = 0.1 is achieved 
for L gain sa 0.34 cm. 

While qualitatively the results reported below are valid 
for general periodic modulations R(rf), the quantitative 
analysis is performed for the semi-infinite lattice of the 
form R(rj) = p r sin 2 (2n), with p r being the depth of the 
lattice proportional to the refractive index modulation, 
which is placed at 77 > (notice that the period of the 



structure can be also made 7r/2 by simple renormaliza- 
tion). At 77 < we have a homogeneous medium where 
R(rf) = 0. We suppose that localized gain, whose profile is 
described by the function 7(7/), is realized in the vicinity of 
a near-surface lattice channel. In the simplest case we con- 
sider 7(77) exactly coinciding with R(ri) in the first channel 
of the lattice, i.e. 7(7/) = pi sin 2 (27/) (pi > is the linear 
gain coefficient) for < r\ < ir/2, and 7(77) = otherwise, 
but the situation when the peak position of 7(77) is shifted 
with respect to first maximum of R(tj) by a distance r] s 
will be considered too. 

Localized dissipative solitons are searched in the form 
QiVtO = w (v) exp[ifo£ + i0(r/)] with w(j]) and 9{rf) being 
the amplitude and the real stationary phase of the field, 
and b being the propagation constant. Exponentially lo- 
calized modes in periodic media |18j and conservative sur- 
face modes [3H5] emerge when the propagation constant b 
falls into one of the gaps of spectrum of periodic guiding- 
structure. It turns out that this is also true in the case of 
dissipative surface solitons. To prove this, we rewrite (TTJ 
in terms of the real functions w{rf) and 11(77) = 9^: 
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and consider first the limit 77 — > —00. Since in this limit 
R = 0, one readily finds the explicit asymptotics 

w = A-e^'" - crA 3 _e 3v/25 "/(86) + C(e 5V2S ")(3a) 
v = -aA 2 _e 2V2 ^'"/(2V2b) + 0{e iV2 ~ bn ) (3b) 

where A_ is a real constant, depending on the total energy 
flow U = J_ w 2 drj. The first (trivial) consequence of the 
obtained asymptotics is that localized solutions exist only 
in domains outside the linear spectrum, i.e. at b > 0. Sec- 
ond, important conclusion is that in the asymptotic region 
v 2 w ~ exp(5v / 2677), i.e. decays faster than w 3 . Thus, at 
77 — > —00 the shape of the mode is described by the conser- 
vative properties of the medium, i.e. 9(77, £) exponentially 
approaches the standard stationary NLS soliton. 

Similar, but more sophisticated, analysis can be per- 
formed for the limit 77 — > 00, where the linear lattice is 
present. One still can prove that v 2 w decays faster than 
u> 3 , and thus the leading order for the field amplitude is 
given by the Floquet theorem: w = A + exp(—fir])P n (r]) 
where P n (v) i s an 7r /2 or 7r periodic function, fj, is the 
Floquet exponent which is determined by the detuning of 
the propagation constant from the band-edge towards the 
n-th stop gap (the lower and upper boundary of n-th gap 
will be designated by b~ and 6+, respectively, the num- 
ber of the first finite gap is set to be n = 1, while the 
semi-infinite gap is denoted by n = 0), and A + is the nor- 
malization constant (see e.g. [18] for more details). We 
illustrate the band-gap spectrum in Fig.QJa). 

Thus, for existence of a localized surface mode one 
has to require that the propagation constant lies beyond 
the allowed band of the linear spectrum of the uniform 
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Fig. 1: (a) Band-gap spectrum of the infinite lattice, (b) b 
vs pi in focusing (curve 1) and defocusing (curve 2) media at 
a — 1. The horizontal lines show the limits of the total gaps. 
U vs pi in focusing (c) and defocusing (d) media. Panels (b)- 
(d) correspond to rj s — 0. (e) U vs tj s at a — 0.5. (f) Maximal 
positive (r}f) and negative (r]~) shifts of amplifying domain 
with respect to R(rj) at which surface soliton still exists vs pi 
at q = 0.5. The circles in (b) and (e) correspond to solitons 
shown in Fig. [2ja) and (b), the circles in (c),(d) correspond to 
solitons shown in the right column of Fig. [2] 



medium, on the one hand, and belongs to one of the gaps 
of the periodic structure on the other hand. The respective 
domains will be termed total gaps (to distinguish them 
form the own gaps of the lattice: as it is clear total gaps 
represent a subset of the lattice gaps). For example, for the 
case of p r = 5, studied below in details, there are only two 
total gaps: the semi-infinite total gap b £ (2.875, oo) and 
the total finite gap (0, 1.840) while the lattice gaps, shown 
in Fig. [Ha), are given by (2.875, oo), (-0.645, 1.840), etc. 
(i.e. 677 = 2.875, 6+ = 1.840, 67/ = -0.645 etc.) 

The established constrains on 6, naturally impose lim- 
itations on possibility of the excitation of the surface 
modes. The mismatch between the boundaries of the gaps 
in the left and right hand structures implies the existence 
of a threshold value U cu t, of the energy flow. In order 
to show this, let us consider a mode whose propagation 
constant b belongs to the total semi-infinite gap and ap- 
proaches 67, i.e. < ii - Jq C 1. Considering r\ > 0, 



we observe that when b — > 67 the amplitude w tends to 
zero: u> max — > (see e.g. [T^] and references therein). On 
the other hand, considering ([2]) at 77 < as an ODE defin- 
ing the shape of the soliton, the smallness of w means that 
the terms v 2 w in ([2"a| can be neglected and the asymptotic 
behavior is described by the conservative NLS equation, 
i.e. by (|5aj) . Since, b does not go to zero (due to finite 
value of 677 ), also the nonlinear term w 3 can be neglected 
and the field behavior is described by the linear equation 
Wrfrf = 2bw. Thus the function w must be exponentially 
decaying, what for the linear ODE at hand is only pos- 
sible if w v (0) = >/2bw(p). As it is clear this is an extra 
condition, in addition to the continuity of q and q^ , which 
must be satisfied at the boundary, i.e. at 77 = 0. In a 
general situation this is impossible with only two avail- 
able constants A + and A_, which are determined by the 
total energy U and by the properties of the linear lattice 
at 77 > 0. In other words, by assuming that the amplitude 
of the mode can go to zero we have arrived at a contradic- 
tion. Thus, there exists a minimal threshold value of ui max 
above which surface modes can exist. Taking now into ac- 
count the Sobolev inequality |iy ma x| 2 < 2U J Iw^d-q we 
conclude that there exists a threshold for the energy flow. 
Moreover, since the limit b — > 6(7 would imply w max — > 
we finally conclude that there exits also a cut-off value 
bcut > b$ such that only for b > b cu t one can find dissipa- 
tive surface modes. Below we illustrate these properties in 
numerical examples. We notice that the above arguments 
are also valid for the pure conservative case, and thus ex- 
plain the threshold values for the energy flow U, observed 
in earlier studies [3HS] . 

Dissipative surface solitons exist not only due to balance 
between diffraction, refraction, and nonlinearity, but also 
due to balance between localized gain and nonlinear losses, 
expressed by the condition a J_ w 4 dr] = J_ ~f{r))w 2 dr]. 
Hence the propagation constant b and the energy flow U 
are determined by the gain pi and by nonlinear losses a. 
Typical dependencies b(pi) for surface solitons in focusing 
medium arc shown in Fig. QI,b). The propagation con- 
stant b in the focusing medium falls into the total semi- 
infinite gap, i.e. 6 > 6(7, and monotonically grows with 
Pi. Such surface solitons exist above some threshold value 
of b, b cut > 67, respectively above the minimal value of 
the gain coefficient, denoted below as p l ° w , and above the 
threshold energy flow U cu t- The energy flow monotoni- 
cally increases with pi everywhere except for a very narrow 
region close to p l ° w [Fig. QJc)]. Typical profiles of dissi- 
pative surface solitons in a focusing medium are shown 
in Fig. HJa) when gain is realized in near-surface lattice 
channel. For low values of pi the surface solitons expand 
considerably into the lattice region and acquire shape rem- 
iniscent to oscillations of shape of Bloch state bordering 
the respective gap edge. Due to energy flow in the trans- 
verse direction and despite the presence of nonlinear losses 
in the entire medium solitons may extend far beyond the 
amplifying region. With increase of gain pi the light grad- 
ually concentrates in the near-surface channel. It should 



p-3 



Y. V. Kartashov, V. V. Konotop, and V. A. Vysloukh 



be stressed that solitons may form in near-surface lattice 
channel even when the gain is displaced by a distance rj s 
with respect to the first maximum of the lattice. Typical 
profile of a surface soliton supported by such shifted gain 
landscape is shown in Fig.^b) - here the maximum of the 
field remains in the surface channel despite the fact that 
gain is realized almost between first and second channels. 
The energy flow of surface solitons first decreases with r] s , 
and then increases when shift approaches maximal value 
beyond which surface soliton can not form in the first lat- 
tice channel [Fig.QJe)]. The maximal possible shift of gain 
landscape quickly increases with pi and saturates already 
at pi = 3 [Fig. QJf)]. Remarkably, surface solitons may 
form not only when gain profile is shifted into the depth of 
the lattice (positive r] 8 ), but also when gain is shifted into 
uniform medium (negative rj s ). When shift rj s becomes 
sufficiently large solitons may form in second, third, etc, 
channels of the lattice. Representative examples of profiles 
of dissipativc surface solitons in second lattice channel are 
shown in Fig. [2fc). Such solitons feature smaller thresh- 
olds (both in terms of pi and U) for their existence than 
their counterparts in the first lattice channel. 

In the case of defocusing medium localized gain can sup- 
port dissipative gap solitons featuring characteristic os- 
cillating tails (inside the lattice) at the surface of semi- 
infinite lattice. The propagation constant of such solitons 
falls into the finite total gap, b £ (0, bf) [see Fig. Q] (b) 
where there is only one total gap] and decreases with p^. 
Like their counterparts in focusing medium, now the soli- 
tons emerging from the finite total gap exist above the 
minimal value of gain coefficient p l ° w (for this value of 
linear gain the propagation constant approaches a cut-off 
value that is close to the upper edge of the total finite gap). 
However now, due to finiteness of the gap, there exists also 
the upper limit for the linear gain p^ pp at which b reaches 
zero value, i.e. dissipative gap solitons can be found for 
plow < p < p" pp . Respectively, the energy flow takes on 
the values from the finite interval, where it is the increas- 
ing function of pi [Fig.QJd)]. When pi — > p l ° w gap surface 
solitons expand dramatically into lattice region, but re- 
main well localized inside uniform medium [Fig. EJa)] , in 
accordance with the asmyptotcs (j3|). The best overall lo- 
calization is achieved for intermediate pi values [Fig. Etc)]. 
When b approaches zero (respectively pi approaches p" pp ) 
the gap surface solitons again become poorly localized 
due to appearance of long tails in the uniform medium 
[Fig. EJf)]. The domains of existence of dissipative sur- 
face solitons in both focusing and defocusing media on 
the (a,pi)-plane are shown in Fig. [3] Minimal gain p l ° w 
required for the existence of solitons in uniform medium, 
as well as the width p™ pp — p l ° w of the band of gain coeffi- 
cients where solitons exist in defocusing medium, increase 
with nonlinear losses a. 

Finally, we analyzed stability of obtained soliton solu- 
tions. We have performed both, the linear stability analy- 
sis and direct propagation method in the presence of input 
perturbations. Both of them showed that the dissipative 
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Fig. 2: Surface solitons in focusing (left column) and defocus- 
ing (right column) media. Panel (a) shows solitons residing in 
first lattice channel at a = 0.4, rj s — 0. Panel (b) shows solitons 
residing in first lattice channel for various shifts of amplifying 
domain with respect to conservative lattice at pi — l,a = 0.5. 
Panel (c) shows solitons residing in second lattice channel at 
a — 0.4, rj s — 7r/2. Gray regions indicate guiding lattice chan- 
nels, while cyan regions show amplifying domains. In (b) we 
show only amplifying domain corresponding to rj a = 0.7. 



surface solitons are exceptionally robust and can with- 
stand even strong shape deformations almost in the en- 
tire existence domain (we considered the lowest branches 
only). 

To conclude, we have demonstrated that an inter- 
face with a gain between a periodic and a homogeneous 
nrcdium can support a diversity of the surface solitons. 
The propagation constants of the modes belong to one of 
the total gaps of the structure. Such modes are character- 
ized by the presence of threshold values of the energy flow 
and cut-off values of the prorogation constant and respec- 
tively of the linear gain. The dissipativc surface solitons 
are attractors and therefore can be easily excited from a 
wide range of initial conditions. Our analysis clearly indi- 
cates that similar modes can be obtained in a more general 
situation where an interface with gain separates two differ- 
ent periodic media. In this last case one can expect larger 
diversity of the surface modes, e especially when the com- 
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Fig. 3: Domains of existence of surface solitons in focusing 
(a) and defocusing (b) media. In the focusing medium surface 
solitons exist for all pi exceeding lower value p l ° w , while in 
defocusing medium surface solitons exist between lower p\ ow 
and upper p" pp values of gain coefficient. 



posite structure is characterized by more that one total 
finite gap. Finally the reported surface dissipative soli- 
tons seem to be very promising objects for spectroscopy, 
sensors, excitations of nano-particles, etc. since their tech- 
nological manufacturing is already available. 
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